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We study the non-equilibrium coherent spin mixing dynamics in ferromagnetic spin-1 and antifer¬ 
romagnetic spin-2 thermal gases of ultracold atoms. Long lasting spin population oscillations 

with magnetic field dependent resonances are observed in both cases. Our observations are well 
reproduced by Boltzmann equations of the Wigner distribution function. Compared to the equation 
of motion of spinor Bose-Einstein condensates, the only difference here is a factor of two increase 
in the spin-dependent interaction, which is confirmed directly in the spin-2 case by measuring the 
relation between the oscillation amplitude and the sample’s density. 
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I. INTRODUCTION 

Ultracold spinor quantum gas has been a subject of 
growing interest in recent years due to its internal degrees 
of freedom which connect it naturally to quantum mag¬ 
netism and many other important physics problems. Ex¬ 
tensive experimental and theoretical investigations have 
been carried out on both the spinor gas itself and its vari¬ 
ous applications, such as ground state phase diagram, do¬ 
main formation, topological excitations, spin squeezing, 
magnetometry and spin mixing dynamics [2 |2]. Cen¬ 
tral to the understanding of these intriguing phenom¬ 
ena is the spin-dependent interaction and its competition 
with other relevant energy scales [3H5]. Non-equilibrium 
coherent spin mixing dynamics, manifested as coherent 
spin population oscillations, is a direct paradigm of the 
spin-dependent interaction and its interplay with Zee- 
man energies [211]. Such dynamics have been observed 
in spin-1 [5HTU] and spin-2 ®^Rb [TT] . and spin-1 ^^Na 
Bose-Einstein condensates (BECs) [12] early on. 

In analyzing the coherent spin mixing dynamics in 
spinor BECs, nice agreements can typically be obtained 
under the single-mode approximation (SMA), in which 
the external spatial and the internal spin degrees of free¬ 
dom are separated from each other [s] 12 [m [n] . It is 
thus natural to ask whether such dynamics can be ob¬ 
served with multi-spatial-mode spinor gases. Although 
collisions in these gases are typically considered random 
or incoherent in the spatial degrees of freedom, coher¬ 
ence among spin degrees of freedom can persist for a 
long time |I3) . Indeed, coherent spin mixing dynamics 
were observed recently in a thermal spin-I ^^Na gas [14] 
and a quantum degenerate Fermi gas of with large 
spin [m [IS]- Results from these experiments were ex¬ 
plained well by dynamics in the spin degrees of freedom 
only, indicating spin and spatial modes can still be largely 
decoupled under right conditions even without conden¬ 
sates. 

In this work, we present experimental investigations of 
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coherent spin mixing dynamics in ultracold thermal ®^Rb 
spinor gases. Different from the antiferromagnetic spin- 
1 Na spinor gas the spin-dependent interaction in 
spin-1 ®^Rb is ferromagnetic and typically much smaller 
in magnitude [THUS]. Amazingly, we still observe robust 
and long-time coherent spin mixing dynamics driven by 
the 70 pK spin-dependent interaction in thermal sam¬ 
ples with a typical temperature of 400 nK. In addition, 
we also observe clean coherent spin mixing dynamics in 
spin-2 ®^Rb thermal gas, which has more spin mixing 
channels and thus richer dynamics [20]El]. In both cases, 
dependences of the spin dynamics on external magnetic 
fields are studied in detail and are explained well by a 
theoretical model developed under a single-spatial-mode 
approximation. 

The rest of this paper is organized as follows. In sec¬ 
tion II, we derive the collisionless Boltzmann equation for 
modeling the spin dynamics of thermal spin-1 and spin-2 
Bose gases. In section III, we describe the experimental 
setup and data taking procedure. In section IV, experi¬ 
ment results are presented and compared with numerical 
simulations based on the model. We then conclude the 
paper in section V. 

II. THEORY 

While spinor BECs are well described by coupled 
Gross-Pitaevskii equations, pure thermal spinor gases 
can be dealt with the semiclassical Boltzmann transport 
equation with the Wigner function as the distribution 
function. Following this approach, several groups have 
predicted the existing of spin waves and spin mixing oscil¬ 
lations in spin-1 thermal gases [22] [23]. These theoretical 
results were successfully applied to the spin-1 Na thermal 
spinor gas in reference M- It was found that under the 
right experimental conditions, the spin dynamics can be 
separated from the multi-spatial modes. Compared with 
spin-1 spinor BECs, one needs only to multiply a fac¬ 
tor of two in the spin-dependent interaction coefficient 
to account for the fact that thermal atoms are distin¬ 
guishable [14]. Here we follow the same formalism but 
generalize it to include both spin-1 and spin-2 cases. 
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In the second-quantization language, the Hamiltonian 
for spin-1 and spin-2 atoms of mass m in an external 
magnetic field B can be expressed jointly as 


action potential U, 

Ukj{r) = [qF^ + 5oTr(n) -h 


H = J dr i)l + Vir) + qF^"^ 


kj 


kjlm 




kj 


where q = quadratic Zeeman energy, which 

is positive for spin-1 and negative for spin-2 ®^Rb atoms. 
Here g is the hyperfine Lande g-factor, is the Bohr 
magneton, and lS.Eh / is the ground-state hyperfine split¬ 
ting. The linear Zeeman energy is gauged out due to 
total spin conservation. V (r) is the external potential, 
F is the vector spin operator and {'4’\) is the bosonic 
field annihilation (creation) operator for spin projection 
k. h is the reduced Planck’s constant. 

The interaction coefficients go, gi, and g 2 can be ex¬ 
pressed in terms of the s-wave scattering lengths ag, 02 , 
and 04 in the total spin 0, 2 and 4 channels of two col¬ 
liding atoms. For the spin-1 case, the total spin can 
be 0 and 2, and only go = 47rIi^(ao -I- 202 )/3m which 
represents the density-dependent interaction, and gi = 
47r/i^(o2 — ao)/3m which represents the spin-dependent 
interaction, are present. The interaction between spin-1 
®7Rb is ferromagnetic with gi negative. For the spin- 
2 case, the total spin can be 0, 2 and 4. Besides 
go = 47rh^(4a2 -I- 3a^/7m and gi = 4Trh?{a4 — a2)/7m, 
there is also a third term ^'2 = 47rIi^(7ao —10a2-l-3a4)/7TO, 
which represents the singlet-pairing interaction. Spin-2 
®^Rb is believed to be antiferromagnetic with gi positive 
and 52 negative but of much smaller magnitude |24j . 

Following the derivation of spin-1 Boltzmann equa¬ 
tion [H121 US], we introduce the Wigner distribution 
function which has matrix elements 

fkj{r,p,t)= /'dr'e"*P'''/^(x/>](r-r 72 ,t) 

J (1) 

7fe(r-f r72,f)). 

From the Wigner function, we can obtain experimen¬ 
tal observables such as the local density nkj{r,t) = 
f dp fkj (r, P, t) / {2TTh)^. Taking the Hatree-Fock approx¬ 
imation, we have the standard collisionless Boltzmann 
equation of the Wigner function 

d 1 

jrJkj + - ■ Vrfkj - VrF • Vpfkj - ^[U, f]kj 

Ub TTl %th ( 2 ) 

-^{VrC/,Vp/}fc, =0, 

where we have introduced matrix elements of the inter- 
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2g2 


Tr (F^n) F^ + giY 


kj 


(3) 




(-) 


'^—j, — k • 


Here g = x, y and z, and Tr is the trace operation. The 
factor of two in the last term is a result of equal Hartree 
and Fock term contributions. Note that this term van¬ 
ishes for spin-1 case. 

In the limit of strong trapping potential where the 
spatial motion is moderately faster than the spin dy¬ 
namics characterized by interaction energies gi_ 2 Ti'(n), 
the spatial dependent interactions can be averaged out. 
Thus we can decouple the spinor evolution cr^,j(t) = 
from the spatial and momentum distri¬ 
bution w{v, p) and express the Wigner function matrix 
elements as fkj{r,p,t) = w{Y,p)akj(t), in analogous to 
the popular SMA in dealing with spin dynamics in BECs. 
Substituting this into Eq. (|^ and integrating over posi¬ 
tion and momentum, we obtain the following equation of 
motion for the coherent spinor dynamics in a thermal gas 

. (-) 


where the spin-dependent interaction potential 


= 


qF^ +gin YTr{F,a)F,+ 




J kj 


- {-f~^ 

' 2g2n g (7—j^—ji,. 


(5) 


Here, n = f (ir[Tr(u(r))]^/7V with N the total atom num¬ 
ber. 

For spin-1, a is represented by 3x3 matrices. With 
the identity = TrA + A + (Ff_,A) Ff_, 

for any 3x3 matrix A, it can be easily shown that 
the two gi terms have equal contributions in Eq. Q. 
Thus the spin-dependent interaction can be summed as 
2gih^ Tr {F^a) F^, which is a factor of two larger than 

k- 

that in spin-1 pure BECs El- 

In the case of spin-2, a is represented by 5x5 matrices 
and the above identity is not true in general. However, 
if cr is constructed from pure state spinor wavefunctions, 
which is the case in our experiment, contributions from 
the two gi terms are again the same and thus the factor of 
two still holds. Combined with the factor of two in the 
g 2 term, the overall spin-dependent interaction is also 
doubled compared with spin-2 pure BECs. 
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FIG. 1. (color online) Coherent spin population oscillation of spin-1 normal gas and its dependence on magnetic field, 

(a) and (b) are exemplary temporal evolution of po at magnetic fields of 0.1 G and 0.18 G, respectively. Red solid curves are 
for eye guiding, (c) and (d) show the magnetic field dependence of the oscillation period (o) and amplitude (o) with a resonant 
feature located near 0.14 G. Error bars here are from fitting of the oscillations. Blue solid curves are fittings with Eq. 


III. EXPERIMENTS 

Our single vacuum chamber experimental setup has 
been descried before [251 ES] . In brief, we prepare the 
ultracold ®^Rb sample by evaporative cooling in a hybrid 
magnetic quadrapole plus optical dipole trap(ODT) |5S]- 
128] first. The magnetic trap ensures a hundred percent 
spin polarized sample in the |1,—1) hyperfine Zeeman 
state. The sample is then transferred to a crossed ODT 
in which the final evaporation is performed within a 
weak magnetic field to preserve the atom’s spin state. 
In the same setup, we can produce a pure ®^Rb BEG 
with 2 X 10® atoms. For the current experiment, we 
control the atom number and stop the evaporation be¬ 
fore the BEG phase transition to make pure thermal 
gases. The measured typical final trap frequencies are 
[lOx, u}y,Wz) = 27 r X (190,211,113) Hz. 

For investigating the spin-1 case, after the final evap¬ 
oration, we hold the sample in the ODT for several 
hundred milliseconds to ensure full thermalization be¬ 
fore the magnetic field is set to a desired value along 
the z-axis, with the transverse magnetic fields compen¬ 
sated to less than 3 mG. To initialize the spin dy¬ 
namics, we directly apply a resonant radio frequency 
(rf) 7r/2-pulse to transfer the sample from |1,—1) hy¬ 
perfine state to the fully transversely magnetized ini¬ 
tial state Ci(0) = (1/2, l/-\/2,1/2)^. Here Ci = 

with /0 -i,o,-hi and 

0_i.o,-i-i as fractional populations and phase in each spin 
components respectively. The system’s magnetization is 
defined as — /o_i. 

To investigate spinor dynamics in the spin-2 case, 
we first transfer the atoms to |2, — 2) hyperfine state 
with a microwave rapid adiabatic passage at a H field 
of 1.3 G with a near unity efficiency. The magnetic 


field is then changed adiabatically to a desired value 
before a fully transversely magnetized state C 2 ( 0 ) = 
(1/4,1/2, -\/3/8,1/2,1/4)"^ is prepared with a rf 7r/2- 
pulse. The magnetization is defined as 2p+2 — 2p_2 + 
p+i — p-i in this case. 

After initialization, the system is in a far-from- 
equilibrium state and spin mixing dynamics will start. 
After holding the sample in the trap for a range of 
time for the dynamics to evolve, the ODT is turned 
off abruptly and atoms in different spin states are de¬ 
tected with the time-of-flight Stern-Gerlach absorption 
imaging technique after 12 ms expansion. Our absorp¬ 
tion imaging setup is calibrated with the high saturation 
method |29l|30]. The number of atoms in each spin state, 
Ni, is extracted from the images and the fractional pop¬ 
ulation Pi = Ni/N is then calculated. 

IV. RESULTS AND DISCUSSIONS 
A. Spin-1 case 

In Fig.Wa) and (b) we present spin population oscilla¬ 
tions for ^Rb thermal spin-1 gases at magnetic fields of 
0.1 G and 0.18 G, respectively. We can immediately see 
that the oscillation depends strongly on magnetic field. 
This can be understood from which has only the 2gi 
term and the quadratic Zeeman q term for spin-1. While 
the negative 2gi term favors the ferromagnetic state, the 
positive q favors the polar state |3]. Non-equilibrium spin 
mixing dynamics is just the result of the competition be¬ 
tween them. Similar to spin-1 BEG, depending on their 
relative strengths, spin dynamics can be divided into the 
interaction regime and the Zeeman regime |S]. With the 
2gi term fixed by the density, external magnetic fields 
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can be applied to tune the system into either regimes. 

When the magnetic field is low and thus q is small, the 
system is in the interaction regime and the oscillation 
period is predominately determined by the 2gi term [5]. 
This is the case in Fig. Sa). At higher magnetic field B, 
q (^725^ Hz/G^) dominates the rather small 2gi term 
(251 (n) ~ 1.45 Hz for our typical density), and the system 
enters the Zeeman regime in which the oscillation period 
is oc 1/g [5]. In both cases, these oscillations last for a 
rather long time, although only data in the first 3 seconds 
are shown. For longer time, the oscillation continues but 
becomes non-periodic. To extract the oscillation period 
and amplitude, we fit the first several oscillations to a 
damped sinusoidal function. 

For the current initial state, the crossover between 
these two regimes happens when q « 2gi (n) mui]. A 
resonance feature is observed near this crossover in our 
experiment, as illustrated in Fig. [^c) and (d). The res¬ 
onance happens at about 0.14 G, evident by the longest 
period and the largest amplitude at this field. On the 
higher field side where the dynamics is dominated by 
q, the oscillations become faster with smaller amplitude. 
This is similar to the detuned Rabi oscillations in driven 
two level systems. Eventually, when the magnetic field is 
too large, which corresponds to the large detuning case, 
the oscillation amplitude becomes too small to be ob¬ 
served. On the lower field side, the behavior is quite 
different. While the amplitude also keeps decreasing, the 
period levels off to ~ l/2gi{n). 

As illustrated by the solid curves in Fig. [^c) and (d), 
these behaviors are well captured by the simulation with 
Eq. Q. These curves are fits to the data points with 
the measured mean density of (n) « 3.0 x lO^^cm”^ 
and a residual magnetization of 0.06(2) due to imperfect 
control of the rf pulse area in the initial state prepara¬ 
tion. With gi as the only fitting parameter, we obtain 
02 — Oq ~ —1.00 ± 0.12 qb {dB is the Bohr radius), con¬ 
sistent with the reported value of 02 — oq « —1.07 as 
in reference m- The rather small but non-zero mag¬ 
netization also explains the non-diverging on resonance 
oscillation period [7]. 

To our knowledge, the current work is the first obser¬ 
vation of magnetically tuned spin oscillation resonance in 
the spin-1 ®^Rb spinor gas. Previous works with BECs 
were performed with either an initial state with large 
magnetization [5] or in a quasi-one-dimensional trap [5] . 
In the former case, the resonance does not exist [7]. For 
the latter, the spin healing length ^ = %/^J2mn\gi \ is 
smaller than the size in the elongated direction and thus 
SMA is violated. Spin mixing dynamics is unstable |32j in 
this case, as perturbations can cause irreversibly conver¬ 
sion of the spin-dependent energy to collective excitation 
modes. This will lead to the formation of multiple spin 
domains and destroy the spin coherence within a single 
full oscillation, making it impossible to observe the reso¬ 
nance. This problem is mitigated in the current work by 
the tight and near 3-D crossed trap. 


B. Spin-2 case 

The interaction between spin-2 ®^Rb atoms is antifer¬ 
romagnetic with gi > 0, 52 < 0, and I 52 I "C gi [24]. Un¬ 
like the spin-1 case, spin-2 spinor oscillations can have 
more than one interaction channel and the spin-2 equa¬ 
tion of motion has no exact analytic solutions. The spin-2 
®’^Rb gas is also subject to inelastic hyperfine changing 
collisions which greatly limit its lifetime to be about 0.5 
s. 

As shown in Fig. [^a) and (b), several full oscillations 
show up within hundreds of milliseconds. The observed 
behaviors are similar with those studied previously on 
spin-2 ®^Rb BEGs with the same initial state mi^ 
and can be understood intuitively from the competition 
between q and the spin-dependent interactions. Indeed, 
ignoring the small 52 term, for the chosen initial state 
C 2 , approximate solutions have been obtained in refer¬ 
ences [in ESI both in the deep interaction and Zeeman 
regimes. These solutions can be directly generalized to 
our case by replacing the gi factor with 2gi , as discussed 
in Section II. 

The factor of two thermal enhancement of the spin- 
dependent interaction can be verified experimentally. In 
the Zeeman regime, spin mixing process (0) -I- (0) 0 
(—1) -I- (1) dominates and the oscillation follows the fun¬ 
damental period ir/q [TTJESj, as depicted by the dashed 
curve in Fig. [^c). Although the amplitude is typically 
small, the oscillations are quite fast so that several pe¬ 
riods can be observed clearly. As shown in the inset of 
Fig. I^c), with the density varied from 1.3x10^^ cm ^ 
to 2.9x10^^ cm“^ at 0.6 G, the periods are nearly con¬ 
stant, while the amplitudes increase linearly with a slope 
of 0.022(2.5)/10^^cm“^. With the best known value 
04 — 02 « 6.95 gb [31], this slop gives the oscillation am¬ 
plitude as 2.0(3) X 3^1 (n) /8q. While for F = 2 spinor 
BEG, the oscillation amplitude is 3gi {n) /8q [TTJ 155] . 
This measurement thus confirms our theoretical predic¬ 
tion that, given a pure initial spin state, the factor of two 
enhancement in the gi term holds for the spin-2 thermal 
gas. 

Fig. [^a) is taken in the interaction regime at 0.1 G, 
where q is smaller than the gi term; while Fig. [^b) is 
taken at 0.45 G with q much larger than the gi term. 
Besides the rather different oscillation amplitudes and 
periods, we also find that in the interaction regime, pop¬ 
ulations only oscillate between mp = 0 and mp = ±2 
states, while those of the mp = ±1 states are nearly 
0.25. In the other regime, all spin states are involved, 
but the mp = ±2 states have a smaller oscillation am¬ 
plitude. In the interaction regime at very low field, the 
oscillation period is 7r/4gi (n) [TTJ|33], which is typically 
long as gi (n) is small. As a result, few oscillations can be 
observed within the sample lifetime and these oscillations 
also show strong damping. 

We have carried out similar measurements with mag¬ 
netic fields ranging from 0.1 G to 0.6 G and extracted 
the amplitudes and periods by fitting the oscillations 
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FIG. 2. (Color online) Spinor dynamics and its dependence on magnetic field for the ®^Rb spin-2 normal gas. (a) and (b) show 
evolutions of po(A), pi(»), p-i(V), p 2 (B), and p_ 2 (<l) at 0.1 G and 0.45 G, respectively. Red solid curves in (a) and (b) are 
fits to the damped sinusoidal function for po, pi, and p 2 only. The measured periods (o) and amplitudes (o) of po vs. magnetic 
fields are summarized in (c) and (d). The accompanying error bars are fitting errors. The blue solid curves here are obtained 
from numerical calculations with Eq.Q and the black dashed curve is a plot of n/q. The inset of (c) shows the dependence of 
the oscillation periods (up triangle) and amplitudes (squares) on densities at 0.6 G (see text). 


to the damped sinusoidal function. As summarized in 
Fig. [|;c) and (d), a resonance is observed near 0.3 G. 
Close to resonance, the approximate solution fails. Nu¬ 
merical calculation of the equation of motion in Eq. Q 
is thus necessary to fully describe the magnetic field de¬ 
pendence. This is performed with the gi term obtained 
from the enhancement factor verification above and the 
experimentally measured number densities. As shown by 
the solid curves in Fig. [^c) and (d), without any free 
parameters the results already agree with our measure¬ 
ments very well. 

V. CONCLUSION 

We have observed and analyzed the non-equilibrium 
interaction-driven collective spin mixing dynamics in fer¬ 
romagnetic spin-1 and antiferromagnetic spin-2 gases of 
ultracold but non-condensed ®^Rb atoms. These dynam¬ 
ics and their magnetic field dependence are proved to be 
the same as those found in spinor BECs under SMA, with 
only a factor of two enhancement in the spin-dependent 


interactions. In the spin-1 case, we can observe oscilla¬ 
tions lasting for a very long time limited only by the trap 
lifetime. Although these oscillations already become ir¬ 
regular at 3 s, we nevertheless cannot tell any obvious 
equilibrium state is reached within 10 s. Spin domains 
formation are fully suppressed in both cases, but even for 
the spin-1 case without hyperfine changing losses, damp¬ 
ing is still observed for most oscillations. This maybe 
come from the ignored collisional integral in the Boltz¬ 
mann equation US] as well as technical imperfections, 
such as the residual magnetic field gradient. 
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